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In the quantum Hall effect (QHE) the differential resistivity r^x = r vanishes within a range 
where the Hall resistivity forms a plateau. A microscopic theory is developed, starting with a 
crystal lattice, setting up a BCS-like Hamiltonian in terms of composite bosons, and using statistical 
mechanical method. The main advantage of our bosonic theory is its capability of explaning the 
plateau formation in the Hall resistivity, which is assumed in the composite fermion theories. In 
the QHE under radiation, the resistivity vanishes within a range with no plateau formation. This 
is shown in terms of two-channels model, one channel excited by radiation where the supercurrents 
run and the other (base) channel in which the normal currents run. The transport diamonds (TD) 
and the zero direct current anomaly (ZCA) occur when the resistivity r is measured as a function 
of magnetic field and direct current (DC). The spiral motion of an electron under a magnetic field 
can be decomposed into two, the cyclotron motion with the cyclotron mass m* and the guiding 
center motion with the magnetotransport M* . The quantization of the motion generates magnetic 
oscillations in the density of states. The magnetoconductivity is calculated, using kinetic theory and 
quantum statistical mechanics. The TR and ZCA are shown to be a breakdown of QHE. The integer 
QHE minima are shown to become the Shubnikov-de Haas (SdH) maxima progressively as the DC 
increases. The ZCA at low temperatures (T = 0.253-1.2 K) is temperature-dependent, which is 
caused by the electron-optical-phonon scattering. 

PACS numbers: 73.43.-f, 73.43.Qt, 74.25.Ha 



I. INTRODUCTION 



Recently, Studenikin et al. [\\ discovered transport 
diamonds (TR) and zero-current anomaly (ZCA) in 
In^Gai_^As/InP and GaAs/AlGaAs. The TR and ZCA 
occur when the differential resistance r^x — dVxx/dl 
of a Hall bar sample after the red-light illumination is 
plotted in the plane of the magnetic field and the di- 
rect current (DC). See Fig. |31 which is reproduced from 
Ref. Fig. 2. Diamond-shaped regions developing from 
SdH minima are called transport diamonds and a sharp 
dip in rxx appearing at a narrow horizontal line at zero 
DC is called a zero-current anomaly. The details of the 
experiments and theoretical backgrounds can be found 
in Ref. [l|. The original authors [l| suggested an in- 
terpretation: the breakdown of the quantum Hall effect 
(QHE). We shall show in the present work that this is 
indeed the case based on the composite (c-)boson model, 
the model originally introduced by Zhang, Hansson and 
Kivelson and later developed by Fujita's group Q . In 
the prevalent theories [1], the QHE is discussed in terms 
of the c-fermions. The formation of the Hall resis- 
tivity plateau where the resistivity vanishes, is assumed, 
however. In our c-boson theory, the plateau formation is 



explained from the first principles. In Sec. |TT1 we review 
our microscopic theory of the QHE. The QHE for a sys- 
tem subjected to a radiation is developed and discussed 
in Sec. IIIII The density of Landau states and statisical 
weight in two dimensions (2D) are calculated in Sec. IIVI 
Shubnikov-de Haas (SdH) oscillations 6] in the magne- 
toconductivity and de Haas-van Alphen (dHvA) oscilla- 
tions in the magnetic susceptibility, are jointly called 
magnetic oscillations. They originate in the oscillations 
in statistical weight. Magnetic oscillations are often dis- 
cussed, using a so-called Dingle temperature Q. But this 
is a phenomenological treatment, which must be avoided. 
Following our previous work ^ , we present a microscopic 
theory of the SdH oscillations for a 2D system in Sec. |Vl 
The TR and ZCA are discussed in Sec. IVll 
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II. INTEGER QUANTUM HALL EFFECT 

If a magnetic field B is applied slowly, then the clas- 
sical electron can continuously change from the straight 
line motion at zero field to the curved motion at a finite 
B (magnitude) . Quantum mechanically, the change from 
the momentum state to the Landau state requires a per- 
turbation. We choose for this perturbation the phonon 
exchange attraction between the electron and the fluxon 
(elementary magnetic fiux). Consider an electron with 
a few fluxons. If the magnetic field is applied slowly, 
the energy of the electron does not change but the cy- 
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clotron motion always acts so as to reduce the magnetic 
fields surrounding the electron. Hence the total energy 
of the composite of an electron and fluxon is less than 
the electron energy plus the unperturbed field energy. In 
other words, the composite (c-)particle is stable against 
the break-up, and it is in a bound (negative energy) 
state. The c-particle is simply a dressed electron car- 
rying Q fiuxons. Q — 1,2, ... . Originally, the c-particle 
was introduced as a composite of one electron attached 
with a number of Chern-Simons gauge objects. [3] These 
objects are neither bosons nor fermions, and hence the 
statistics of the composite is not clear. The basic parti- 
cle property (countability) of the fiuxons is known as the 
fiux quantization, see Eq. (|II.15|) . We assume that the 
fiuxon is an elementary fermion with zero mass and zero 
charge, which is supported by the fact that the fiuxon, 
the quantum of the magnetic field B, cannot disappear 
at a sink unlike the bosonic photon, the quantum of the 
electric field E ^]. Fujita and Morabito [10] showed that 
the center-of-mass (CM) of a composite moves follow- 
ing the Ehrenfest-Oppenheimer-Bethe's (EOB) rule [ll| : 
the composite is fermionic (bosonic) if it contains an odd 
(even) number of elementary fermions. Hence the quan- 
tum statistics of the c-particle is established. 

At the Landau level (LL) occupation number, also 
called the filling factor, v — \/Q, Q odd, the c-bosons 
with Q fiuxons are generated and can condense below cer- 
tain critical temperature Tc. The Hall resistivity plateau 
is caused by kind of the Meissner effect as explained later. 

We develop a theory for GaAs/AlGaAs heterojunc- 
tion, the theory which can also be applied to InGaAs/InP 
quantum well. GaAs forms a zincblende lattice. We as- 
sume that the interface is in the plane (001). The Ga'^"'" 
ions form a square lattice with the sides directed in [110] 
and [110]. The "electron" (wave packet) with a negative 
charge — e will then move isotropically with an effective 
mass nil. The As^~ ions also form a square lattice at 
a different height in [001]. The "holes", each having a 
positive charge (-l-e), will move similarly with an effec- 
tive mass 7712. A longitudinal phonon moving in [110] 
or in [110] can generate a charge (current) density varia- 
tion, establishing an interaction between the phonon and 
the electron (fiuxon). If one phonon exchange is consid- 
ered between the electron and the fiuxon, a second-order 
perturbation calculation establishes an effective electron- 
fiuxon interaction 



eliminated: 
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where q (tkoq) is the phonon momentum (energy), Vq{V^) 
the interaction strength between the electron (fiuxon) 
and the phonon. If the energies (£|p_|_q|, Cp) of the final 
and initial electron states are equal, the effective interac- 
tion is negative (attractive) as seen from Eq. (III. II) . 

Following Bardeen-Cooper-Schrieffer (BCS) [l3], we 
start with a Hamiltonian H with the phonon variables 
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where n^^ is the number operator for the "electron" (1) 
hole" (2), fiuxon (3)] at momentum k and spin s with 

We represent the "electron" ("hole") 

'^ks'^'^ks ' where c{c') are annihilation (cre- 
ation) operators satisfying the Fermi anticommutation 
rules: 
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the energy e^-'^^ 
number n'^^ by c'^}^ c'^} , 



Jj)t„(«) 
Vs''^ks 



= 0. 



(H.3) 
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We represent the fluxon number n^^^ by aj^^aks, with 
a(a^), satisfying the anticommutation rules, (|II.3p . 
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The prime on the summation in Eq. pi.2p means the 

restriction: < e):^ < huju, lou = Debye frequency. 
If the fiuxons are replaced by the conduction elec- 
trons ( "electrons" , "holes" ) our Hamiltonian H is re- 
duced to the original BGS Hamiltonian, Eq. (24) of 
Ref. [T3]. The "electron" and "hole" are generated, 
depending on the energy contour curvature sign p^ . 
For example, only "electrons" ("holes"), are generated 
for a circular Fermi surface with the negative (posi- 
tive) curvature whose inside (outside) is filled with elec- 
trons. Since the phonon has no charge, the phonon 
exchange cannot change the net charge. The pair- 
ing interaction terms in Eq. (III. 21) conserve the charge. 

The term -I'oS^'q.S^q',, where vo = \VqV^\ {TiiOoAr^ 
sample area, is the pairing strength, generates 
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the transition in the "electron" states. Similarly, the 
exchange of a phonon generates a transition in the 

(2) (2)i 

"hole" states, represented by — ''^o^k'qs^kq* ' "^^^ 
phonon exchange can also pair-create and pair-annihilate 
"electron" ("hole")-fluxon composites, represented by 

^kqs: -^oB^'is^l^l- At K, the system 
can have equal numbers of — (-I-) c-bosons, "electrons" 
("hole") composites, generated by "I'o-Bk'qs^kqs ■ 

The c-bosons, each with one fluxon, will be called the 
fundamental (f) c-bosons. At a finite temperature, there 
are moving (non-condensed) fc bosons. Their energies 
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are obtained from p^ : 



u;(^")vl/(k,q) =£[^V^|Vl/(k,q) ~ {27Th)-'v* 

X y" d2fc'^(k'^q)^ (n.5) 

where ^'(k, q) is the reduced wavefunction for the fc- 
boson; we neglected the fluxon energy. The Vq denotes 
the strength after the ladder diagram binding, see below. 
For small q, we obtain a solution of Eq. (III.5[) as 



w^-" = Wo + (2/7r)t;^^^q, Wq 



exp(i;*I?o)-i - 1 
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where Wp = {2eF/mj)^^'^ is the Fermi velocity and Dq = 
D{eF) the density of states per spin. The brief derivation 
of Eqs. pi.5l) and (|II.6p is given in Appendix [X] Note 
that the energy Wq depends linearly on the momentum 
magnitude q. 

The system of free fc-bosons undergoes a Bose-Einstein 
condensation (BEC) in 2D at the critical temperature 



k^Tc — 1.24?iWF«y^. 



(11.7) 



A brief derivation of Eq. (III.7|) is given in Appendix |B] 

1/2 

The interboson distance Rq = Uq calculated from this 
expression is 1.247iwF(fcBTc)~^- The boson size rg cal- 
culated from Eq. (jll.6[) . using the uncertainty relation 
(gmaxT-o ~ ^) and l^ol ~ ^bTc is (2/7r);iwF(fcBTc)~\ 
which is a few times smaller than Rq. Hence, the fc- 
bosons do not overlap in space, and the model of free 
bosons is justified. For GaAs/AlGaAs, m* = 0.067 me, 
nie = electron mass. For the 2D electron density 
lO^^cm"^, we have vf — 1.36 x lO^cms"^. Not all 
electrons are bound with fluxons since the simultaneous 
generations of ± fc-bosons is required. The minority car- 
rier ("hole") density controlls the fc-boson density. For 
uq — 10^" cm^^, Tc = 1.29 K, which is reasonable. 

In the presence of Bose condensate below Tc the un- 
fluxed electron carries the energy 15 1 



E 



.(/)2 
k 



(II.8) 

where the quasi-elcctron energy gap A is the solution of 
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Note that the gap A depends on the temperature T. At 
the critical temperature there is no Bose condensate 
and hence A vanishes. 

Now the moving fc-boson below Tc has the energy Wq 
obtained from 

zi«vl/(k,q) =i?f^^q|Vl/(k,q) - {27rhr^v* 



J d2fc'$(k',q). (11.10) 



We obtain after solving Eq. (Ill.lOP : 
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= wo + {2/7:)vl^\ = wo+eg + {2/7r)v'^>q, (11.11) 
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where Wo(T) is determined from 



1 = Dqi^o 



de 



KI + (e^ + A2)V2 



The energy difference: 



Wo{T) - Wo = £g{T) 



(11.12) 
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represents the T-dependent energy gap. The energy Wq 
is negative. Otherwise, the fc-boson should break up. 
This limits eg{T) to be \wo\ at K. The fc-boson energy 
gap £g declines to zero as the temperature approaches Tc 
from below. 

The fc-boson, having the linear dispersion pi. Ill) can 
move in all directions in the plane with the constant 
speed (2/7r)wp^ as seen from Eq. (jll.lll) . The su- 
percurrent is generated by the ± fc-bosons condensed 
monochromatically at the momentum directed along the 
sample length. The supercurrent density (magnitude) J, 
calculated by the rule: (charge e*) x (carrier density tiq) 
X (drift velocity w^), is 



J = e*noVd = e*no(2/7r) 



,(1) 
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(11.14) 



The induced Hall field (magnitude) E'h equals VdB. The 
magnetic flux is quantized 



BA = n^{h/e), = fluxon density. 



Hence, we obtain 
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(11.15) 



(11.16) 



If = riQ valid at = 1, we obtain pu = h/e'^ in 
agreement with the plateau value observed. 

The model can be extended to the integer QHE at 
i'==P,P = l,2,.... The field magnitude is less. The LL 
degeneracy eBA/h is linear in B, and hence the lowest P 
LL's must be considered. The fc-boson density no per LL 
is the electron density Ue over P and the fluxon density 
is the boson density uq over P: 



no = rie/P, n^ = no/P. 



(11.17) 



At = 1/2 there are c-fermions, each with two fluxons. 
The c-fermions have a Fermi energy. The ± c-ferniions 
have effective masses. The Hall resistivity pn has a B- 
linear behavior while the resistivity p is finite. In our 
theory the integer P is the number of the LL's occupied 
by the c-fermions. 

Our Hamiltonian in Eq. pi.2p can generate and stabi- 
lize the c-particles with an arbitrary number of fluxons. 
For example, a c-fermion with two fluxons is generated 
by two sets of the ladder diagram bindings, each between 
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the electron and the fluxon. The ladder diagram binding 
arises as follows. Consider a hydrogen atom. The Hamil- 
tonian contains kinetic energies of the electron and the 
proton and the attractive Coulomb interaction. If we re- 
gard the Coulomb interaction as a perturbation and use 
a perturbation theory, we can represent the interaction 
process by an infinite set of ladder diagrams, each lad- 
der step connecting the electron line and the proton line. 
The energy eigenvalues of this system is not obtained by 
using the perturbation theory but they are obtained by 
solving the Schrodinger equation directly. This example 
indicates that a two-body bound state is represented by 
an infinite set of ladder diagrams and that the binding 
energy (the negative of the ground-state energy) is cal- 
culated by a non-perturbative method. 

Jain introduced the effective magnetic field Q 

B* = B - = B ~ {l/v)neih/e) (11.18) 

relative to the standard field for the composite (c- 
)fermion. We extend this idea to the bosonic (odd- 
denominator) fraction. This means that the c-particle 
moves field-free at the exact fraction. The movement of 
the guiding centers (the CM of the c-particle) can occur 
as if they are subjected to no magnetic field at the exact 
fraction. The excess (or deficit) of the magnetic field is 
simply the effective magnetic field B* . The plateau in 
Ph is formed due to kind of the Meissner effect. Con- 
sider the case of zero temperature near v = 1. Only the 
system energy E matters. The fc-bosons are condensed 
with the ground-state energy wq, and hence the system 
energy E at ly — 1 is 2NoWo, where Nq is the number of 
— fc-bosons (or + fc-bosons). The factor 2 arises since 
there are ± fc-bosons. Away from v = I, we must add 
the magnetic field energy {2fio)~^A{B*)^, so that 

E = 2NoWa + {2fioy^A{B*f. (11.19) 

When the field is reduced, the system tries to keep 
the same number A'o by sucking in the flux lines. 
Thus the magnetic field becomes inhomogeneous out- 
side the sample, generating the extra magnetic field en- 
ergy (2/xo)^^A(i?*)^. If the field is raised, the system 
tries to keep the same number Nq by expeling out the 
flux lines. The inhomogeneous fields outside raise the 
field energy by (2/xo)~^A(i?*)2. There is a critical field 
B* = (4/xo|wo|)^^^- Beyond this value, the supercon- 
ducting state is destroyed, which generates a symmetric 
exponential rise in the resistance R. In our discussion 
of the Hall resistivity plateau we used the fact that the 
ground-state energy wq of the fc-boson is negative, that 
is, the c-boson is bound. Only then the critical field 
B* — {4fj,\wo\y^^ can be defined. Here the phonon ex- 
change attraction played an important role. The repul- 
sive Coulomb interaction, which is the departure point of 
the prevalent fermionic theories [J, [Hj , cannot generate a 
bound state. 

In the presence of the supercondensate, the non- 
condensed c-boson has an energy gap Eg. Hence, the 



non-condensed c-boson density has the activation energy 
type exponential temperature-dependence: 

exp[-e<,/(fcBr)]. (H.20) 

Some authors argue that the energy gap Eg for the integer 
QHE is due to the LL separation = Hluq. But the separa- 
tion HlOc is much greater than the observed Eg. Besides, 
from this view one cannot obtain the activation-type en- 
ergy dependence. 

The BEC occurs at each LL, and therefore the c-boson 
density no is smaller for high P, see Eq. (III.17p . and the 
strengths become weaker as P increases. The most signif- 
icant advantage of our bosonic theory is that we are able 
to explain why the plateaus in the Hall resistivity is de- 
veloped when the resistivity is zero as the magnetic field 
is varied. This plateau formation is phenomenologically 
assumed in the fermionic theories. [4, M 



III. QUANTUM HALL EFFECT UNDER 
RADIATION 

The experiments by Mani et al. (l6j and Zodov et al. 
[l7j indicate that the applied radiation excites a large 
number of "holes" in the system. Using these "holes" 
and the preexisting "electrons" the phonon exchange can 
pair-create ± c-bosons, that condense below Tc in the 
excited (upper) channel. The c-bosons condensed with 
the momentum along the sample length are responsible 
for the supercurrent. In the presence of the condensed c- 
bosons, the non-condensed c-bosons have an energy gap 
Eg, and therefore they are absent below Tc- The fermionic 
currents in the base channel cannot be suppressed by 
the supercurrents since the energy levels are different be- 
tween the excited and base channels. These c-fermions 
contribute a small normal current. They are subject to 
the Lorentz force: F = g(E + v X B), and hence they 
generate a Hall field £'h proportional to the field B. This 
is the main feature difference from the usual QHE (under 
no radiation). 

In the neighborhood of the QHE sX v — 1, the current 
carriers in the base and excited channels are, respectively, 
c-fermions and condensed c-bosons. The currents are ad- 
ditive. We write down the total current density j as the 
sum of the fermionic current density j/ and the bosonic 
current density jf,: 

3 = if + ib = en/w/ -I- eubVb, (HI.l) 

where Vf and Vb are the drift velocities of the c-fermions 
and c-bosons. The Hall fields £'h are additive, too. Hence 
we have 

E^ = Enj + E^^b = VfB + VbB. [111.2) 

The Hall effect condition [E-n — VdB) applies separately 
for the c-fcrmions and c-bosons. We therefore obtain 

j nfVf+nbVbe 
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Far away from the midpoint of the zero-resistance stretch, 
the c-bosons are absent and hence the Hall resistivity i?H 
becomes B/{enf): 



i?H = B /{euf) (far away), 



(III.4) 



after the cancellation oivf. At the midpoint the c-bosons 
are dominant. Then, the Hall resistivity i?H is approxi- 
mately equal to h/e^ since 



Eb_ ^ VbB 

j eubVb 



h 



(midpoint), (HI. 5) 



where we used the flux quantization [B = {h/e)n^], and 
the fact that the flux density equals the c-boson den- 
sity Tib aX V — 1. The Hall resistivity i?H — E-n/ j is not 
exactly equal to h/e^ since the c-fermion current den- 
sity enfVf is much smaller than the supercurrent density 
eUbVb, but it does not vanish. In the horizontal stretch 
the system is superconducting, and hence the supercur- 
rent dominates the normal current: eubVb ^ enfVf. The 
deviation Ai?H is, using Eq. piI.Sp . 



Ai?H 



VfB + VbB 



e{nfVf 
UbB ^ 
eUbVb ~ 



' nbVb) 
h 

72- 



VfB_ 
eufVf 



(IH.i 



If the field B is raised (or lowered) a little from the mid- 
point, Ai?H is a constant {h/e^) due to the Meissner 
effect. If the field is raised high enough, the supercon- 
ducting state is destroyed and the normal current sets 
in, generating a finite resistance and a vanishing Ai?H- 
Hence the deviation Ai?H and the diagonal resistance 
Rxx are closely correlated as observed by Mani [l6j . 

In Fig. 2 of Ref. [l6| (not shown here) , we can see that 
in the range, where the SdH oscillations are observed 
for the resistance without radiation, the signature of os- 
cillations also appear for the resistance Rxx with radi- 
ation. The SdH oscillations arise only for the fcrmion 
carriers. This SdH signature in Rxx should remain. Our 
two-channel model is supported here. 

Mani et al. , Fig. 2 of Ref. JJj , shows that the strength 
of the superconducting state does not change much. The 
2D density of states for the conduction electrons asso- 
ciated with the circular Fermi surface is independent of 
the electron energy, and hence the number of the excited 
electrons is roughly independent of the radiaton energy 
(frequency). The "hole" -like excitations are absent with 
no radiation. We suspect that the "hole" -band edge is 
a distance Eq away from the system's Fermi level. This 
means that if the radiation energy huj is less than Eq, 
the radiation can generate no superconducting state. An 
experimental confirmation is highly desirable here. 

If a bias voltage is applied to the system, then a normal 
current runs in the base channel. In the upper channel 
the supercurrent still runs with no potential drop. The 
both currents run in the same sample space. The ap- 
parent discrepancy in the electric potential here may be 



resolved by considering a static charge Q developed in the 
system upon the field application. That is, the system 
will be charged and the static potential 
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(111.7) 



where C is the system's capacitance, can balance the to- 
tal electric potential while the charging does not affect 
the supercurrents. This effect may be checked by experi- 
ments, which is a critical test for our two-channel model. 

In summary, the QHE under radiation is the QHE at 
the upper channel. The condensed c-bosons generate a 
superconducting state with a gap Eg in the c-boson energy 
spectrum. The supercondensate suppresses the c-particle 
currents in the upper channel, but cannot suppress the 
normal currents in the base channel. Thus, there is a fi- 
nite resistive current accompanied by the Hall field. This 
explains the i3-linear Hall resistivity. 



IV. THE DENSITY OF STATES AND 
STATISTICAL WEIGHT 

We calculate magnetic oscillations in the statistical 
weight for a 2D electron system. Let us take a dilute sys- 
tem of electrons moving in the plane. Applying a mag- 
netic field B perpendicular to the plane, each electron 
will be in the Landau states with the energy given by 

E = iNi^ + l/2)huj„ A^L = 0,1,2,.... (IV.l) 

The degeneracy of the Landau level (LL) is 

eBA/2'nfi, A = sample area. (IV.2) 

The weaker the field the more LL's, separated by hujc, are 
occupied by the electrons. The electron in the Landau 
state can be viewed as circulating around the guiding 
center. 

We introduce kinetic momenta 

^x^Px + eAx, Ily=Py + eAy, (IV.3) 
in terms of which the Hamiltonian T-L for the electron is 



n 



1 



1 



(tf + if) = — n 



(IV.4) 



The vector potential A = (1/2)B X r can be written 
as Ax = -{l/2)By, Ay = {l/2)Bx, A^ = 0. Using 
the quantum condition [a;,pa;] = [y^Py] = [^jU] = 
[PkiPj/] — 0, we obtain 

[Ux,Uy] = -{eh/i)B. (IV.5) 

If we introduce 

(m*)i/2n, = P, (eB)-i(m*)i/2n, = Q, (IV.6) 
we obtain 

n = {l/2)[P^+ujlQ^], (IV.7) 
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FIG. 1: A 2D Landau state is represented by the ring (shaded 
area) of the phase-space volume 27rnAn. 



and 



[Q,p] = in. 



(IV.. 



Hence, the energy eigenvalues are given by {Ni^+l/2)huJc, 
confirming Eq. (jlV.ip . After simple calculations, we ob- 
tain 



dxdUxdydlly — dxdpxdydpy. 



(IV.9) 



We can then represent quantum states by small quasi- 
phase space cell of the volume dxdUx dydHy . The Hamil- 
tonian H in Eq. (jIV.4|) does not depend on the position 
{x, y). Assuming large normalization lengths (Li, L2), we 
can represent the Landau states by the concentric shells 
of the phase space having the statistical weight 



27rnAn • LiL2{2nhy 



eBA 
27rh ' 



(IV.IO) 



where A = L1L2 and huc = A(lP/2m *) = HAn/m*. 
Hence, the LL degeneracy is given by Eq. (jIV.2p . Figure[T] 
represents a typical Landau state in the Ha;-Hy space. 

As the field B is raised the separation TiuJc increases, 
and the quantum states are bunched together. As a result 
of the bunching, the density of states J\f{e) should change 
periodically. 

The electrons obey the Fcrmi-Dirac statistics. Consid- 
ering a system of free electrons, we define the Hclmholtz 
free energy J- by 



J'=Nn-2kBT 



l+e 



(IV.ll) 



where /i is the chemical potential and the factor 2 arises 
from the spin degeneracy. The chemical potential /i is 
determined from the condition 



0. 



(IV.12) 



The total magnetic moment M for the system can be 
found from 



M 



'dB' 



(IV.13) 



Equation (jIV.13p is equivalent to the usual condition that 
the total number of the electrons, iV, can be obtained in 
terms of the Fermi distribution function 



7V = 2^/(i?,). 



(IV.14) 



The LL Ei is characterized by the Landau oscillator 
quantum number TVl. Let us introduce the density of 
state 



dW/dE = M{E) 



(IV.15) 



such that Af{E)dE ~ the number of states having an 
energy between E and E + dE. Wc write Eq. (jIV.llI) in 
the form 



r°° dW 

T = Nn- 2kBT / dE— In 
Jo "'E 

POO 



1 + e^^"-^)/''''^^ 

Nfi-2kBT dEW{E)f{E). (IV.16) 
Jo 



The statistical weight (number) W is the total number 
of states having energies less than 

E ^ {Ni^ + l/2)hL0c. (IV.17) 
For a fixed pair {E, the density of states is 

dW = i^27rHAne[i; - (A^L + l/2)?iwc], (IV.18) 

where Q{x) is the Heaviside step function: 

e(.) = {ji[^>2 • (1^-19) 

We sum Eq. pV.lSp with respect to A^l and obtain 

00 

W(E) = C{nujc)2 Q[e - (2Al + l)7r], (IV.20) 

C ^ 2TTm* A{2ttK)-'^ , e = 2TTElhuj^. (IV.21) 
We assume a high Fermi-degeneracy such that 

^ ~ ep > ^iwc. (IV.22) 



The sum in Eq. (jIV.20[) can be computed by using Pois- 
son's summation formula [l8| 

J2 /(27rn) = — ^ / dTf{T)e-^-\ (IV.23) 



m— — 00 



We then obtain [19| 

W{E) = Wo + W'osc, 
= A{m' l-Kr?)E 



(IV.24) 
(IV.25) 



= Chu..- E ^ Sin . (IV.26) 
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The detailed calculations leading to Eqs. (|IV.24|) - (IIV.26|) 
are given in Appendix \C\ Only the first term v ~ 1 in 
Eq. (jIV.26P will be important in practice for weak fields 



£f > ?jWc, 



(IV.27) 



which will be shown later. 

The term Wq, which is independent of B, gives the 
weight equal to that for a free electron system with no 
field. Note that there are no Landau-like diamagnetic 
terms proportional to the squared field B^. 



V. SHUBNIKOV-DE HAAS OSCILLATIONS 

Let us, first, consider the case with no magnetic field. 
We assume a uniform distribution of impurities with den- 
sity nj. We introduce a momentum distribution func- 
tion (j){p,t), defined such that (j){p,t)(Pp gives the rela- 
tive probability of finding an electron in the element (Pp 
at time t. This function will be normalized such that 



(27r?i)2 



d^p(f>ip,t) 



N 
1 



(V.l) 



where the factor 2 is due to the spin degeneracy. 

The electric current density j is given in terms of (/>(p, t) 

as 



J = 



-2e 



J d^p(l>{p,t)p. 



(V.2) 



The distribution function 0(p, t) can be obtained by solv- 
ing the Boltzmann equation for the stationary homoge- 
neous state, dropping t: 

eE. |-0(p) = ^ / dni{p,9mp')-<l>{p)]p, (V.3) 
op m* J 

where 9 is the scattering angle, that is, the angle between 
the initial momentum p and the final one p', and I{p, 9) 
is the differential cross section. Solving Eq. (|V.3[) . we 
obtain the conductivity as 13 1 



{2T:hy m* 



d^p 



E 



d[_ 

dE J T{E) 



E 



2m* 



where F is the energy (_E)-dependent relaxation rate 



(V.4) 



T{E)^ni I dni(p,9)(l-cos9) — . 

J m* 



The Fermi distribution function 

f{E) EE re^(-E-M) + 1 

is normalized such that 



(V.5) 



(V.6) 



/ dEv{E)f{E), v{E) 
Jo 



Am 

A ' 



(V.7) 



where I'iE) is the density of states per area. We can 
rewrite Eq. (|V.4p as 



e- df\ E 



o2 poo 

cr = / dEiy(E) , , , , , . 

m* Jo ^ \ dEj T{E) 



(V.8) 



The Fermi distribution function f{E) drops steeply 
near E = fi at low temperatures: ksT <C Ef (Fermi en- 
ergy). If the density of states varies slowly with energy 
E, then the delta-function replacement formula 



dl_ 
dE 



5{E ~ ^i) 



(V.9) 



can be used. Using 

dEN{E) E^ dEM{E)f{E), (V.IO) 

and comparing Eqs. (jV.8[) and the Drude formula 



a — n — , 

m* 70 



we obtain 



n r°° 1 



(V.ll) 



(V.12) 



Note that the temperature dependence of the relaxation 
rate 7o(T) is introduced through the Fermi distribution 
function f{E). 

Let us now consider a field-dressed electron (guiding 
center). We assume that the dressed electron is a fermion 
with magnetotransport mass M* and charge e. The ki- 
netic energy is represented by 



n' 



2M* ^ ^ y' 2M* 



(V.13) 



We introduce a distribution function (/?(II, t) in the 
HajIIy-space normalized such that 



(2^ /'^^n^Cn^n,,^ 



N 
~A 



(V.14) 



The Boltzmann equation for a homogeneous stationary 
state of the system is 



e(E + v X B) 



dip 

9n 



j dr!^n,/(n,0)[(p(n')-(p(n)]. 



(V.15) 

where 9 is the scattering angle, that is, the angle between 
the initial and final kinetic momenta (11,11'). In the 
actual experimental condition, the magnetic force term 
can be neglected. Assuming this condition, we obtain the 
same Boltzmann equation (jV.Sp for a field-free system 
except the mass difference. Hence, we obtain 



2e' 



M*{2'Khf 



E 



df 



(V.16) 
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As the field B is raised, the separation fiujc becomes 
greater and the quantum states are bunched together. 
The statistical weight W contains an oscillatory part, see 
Eq. (IRTMI) 



27re' 



n 



'2 



2m* 



(V.17) 



Physically, the sinusoidal variations in Eq. (I IV. 261) arise 
as follows. From the Heisenberg uncertainty principle 
(phase space consideration) and the Pauli exclusion prin- 
ciple, the Fermi energy sp remains approximately un- 
changed as the field B varies. The density of states is 
high when ep matches the A^L-th level, while it is small 
when Ep falls between neighboring LLs. 

If the density of states, A/'(e), oscillates violently in 
the drop of the Fermi distribution function /(e) = 
|^g/3(e-M) _|_ one cannot use the delta-function re- 

placement formula (jV.9l) . The width of df /de is of the 
order k^T. The critical temperature Tc below which the 
oscillations can be observed is 



(V.18) 



Below Tc, we may proceed as follows. Let us consider the 
integral 



/2t:E\ „ n2 



E : 



2M* 



(V.19) 



We introduce a new variable C = /3{E — v)^ and extend 
the lower limit to — oo (low temperature limit): 



1 _ 1 1 

^ ^ / dC----^^. (V.20) 



s'^ + 1' 



With the help of the integral formula 

dC- 



-|- 1 i sinh ira ' 



(V.21) 



which is proved in Appendix [Dl we obtain from 
Eq. dVl9l) : 

TTkBTcos{27reF/hwa) ^2) 



Here, we used 



sinh{27T^M*kBT/heB)' 



M*ii{T = 0) = to*£f. 



(V.23) 



since the Fermi momentum is the same for both dressed 
and undressed electrons. For very low fields, the oscilla- 
tion number in the range k^T becomes great, and hence 
the sinusoidal contribution must cancel out. This effect 
is represented by the factor [smh{2'K'^M*kB,T /heB)]~'^ . 



We now calculate the conductivity, starting with 
Eq. (IV. 8p . For the field- free case, we may use Eqs. (|V.9p 
and (|V.10P to obtain 



n v(e-p)s-p 
70 



ro(eF) 



(V.24) 



For a finite B, the non-oscillatory part (background) con- 
tributes a similar amount: 



n i'{e'p)sF 

r(£F) ' 
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(V.25) 



calculated for the dressed electrons. The oscillatory part 
can be calculated by using the integration formula / in 
Eqs. (IV.19P and (|V.22j) . This part is much smaller than 
v{£-p)e-p /y{s-p) in Eq. (jV.25|) . since the contribution is 
limited to the small energy range k-oT. It is also small 
by the sinusoidal cancellation. We, therefore, obtain 



(1 



7 r(£F) 

_ irk^T cos(2ttsf /hujc) 

" ep sinh(27r2M*fcBr/neB) ' 



(V.26) 
(V.27) 



Strictly speaking, the contribution of the terms 
with ly = 2,3, ••• in the sum Wosc in Eq. (|IV.26P 
should be added. But this contribution, which carries 
[smh{2Tr^iyM* ksT /heB)]-^ , is small since 



smh{2TT^ M*kBT/heB) > 1. 



(V.28) 



In the present theory, the two masses m* and AI* are 
introduced naturally corresponding to the two physi- 
cal processes: the cyclotron motion of the electron and 
the guiding center motion of the dressed electron. The 
dressed electron is the same entity as the c-fermion with 
two fluxons in the QHE theory. 

In summary, the magnetoconductivity cr(_B), given by 
Eq. (jV.16p . may be written out as 



(1 + 0). 



M* 7 M* T{sf) 
In contrast, the conductivity ctq zero field is 
2 n i'{eF)eF 



TO* 70 m* ro(eF) 



(V.29) 



(V.30) 



where we have assumed that the Fermi energy ep remains 
the same for both cases. We note that the magnetocon- 
ductivity a does not approach the conductivity ctq h^ the 
low field limit. In fact, we obtain in this limit (0 = 0): 



(T — (To = e 1^ 



1 



1 



Af*7 TO*7o 



(V.31) 



The difference arises from the carrier difference. 

If the "decay" rate 6 = 2Tr^M*kBT/he defined through 



smh{S/B) = smh{27r^ M*kBT/heB 



(V.32) 
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FIG. 2: The resistance Rxx versus the reduced inverse mag- 
netic field, B'~^ . See Mani [l9[ for the actual reduction. The 
number N in the abscissas is the intersection number between 
the curves with (w/) and without (w/o) microwaves. 



is measured carefully, the magnetotransport mass M* 
can be obtained directly through 



M* = enS/{27T^kBT). 



(V.33) 



Mani measured the SdH oscillations in GaAs/AlGaAs 
[lil, Fig. 1, T = 0.7K. His data are reproduced in Fig.H 
Clearly, we see the diagonal resistance Rxx linearly de- 
creasing with in the low field limit. For high purity 
samples at very low temperature (~ 0.7K), the impurity 
and phonon scatterings are negligible. By the energy- 
time uncertainty principle the dressed electron can spend 
a short time at the upper LL and come back to the ground 
LL with a different guiding center, which causes a guid- 
ing center jump. We assume that the relaxation rate 7 
is the natural linewidth arising from the LL separation 
divided by h, that is, the cyclotron frequency ujc' 

j^ujc = eB/m*. (V.34) 

This generates the desired dependence for Rxx- 
We fitted Mani's data in Fig. H with 



Rxx = A + Bx 



[E cos{27tCx) + F]3 
sinh(£'a;) 



(V.35) 



where A = 2.3, B = -0.18, C = 23.0, D = 3.1, E ^ 22.0, 
and F = 7.0. The fits agree with the data within the 
experimental errors. Using Dx = S/B, we obtain 



M* = 0.30 me, 



(V.36) 



where rrie is the gravitational electron mass. If m* = 
0.067 me, then M* /m* = 4.5. These are reasonable num- 
bers. 

The relaxation rate 7 = r(£F) can now be obtained 
through Eq. (|V.11I) with the measured magnetoconduc- 
tivity. All electrons, not just those excited electrons near 
the Fermi surface, are subject to the electric field. Hence, 
the carrier density n appearing in Eq. (jV.29p is the total 



density n of the dressed electrons. This n also appears 
in the Hall resistivity expression 



PH 



3 



VdB 



envd 

where the Hall effect condition: 



B 

en ' 



En = VdB, Vd = drift velocity 



(V.37) 



(V.38) 



was used. 

The dressed electrons are there whether the system is 
probed in equilibrium or in nonequilibrium as long as 
the system is subjected to a magnetic field. Hence their 
presence can be checked by measuring the susceptibility 
or the heat capacity of the system. All (dressed) elec- 
trons are subject to the magnetic field, and hence the 
magnetic susceptibility x is proportional to the carrier 
density n although the x depends critically on the Fermi 
surface. We shall briefiy discuss the magnetic moment 
and susceptibility. 

The magnetization A^, that is, the total magnetic mo- 
ment per unit area, can be obtained from 



M 



dB' 



(V.39) 



Using Eqs. (|IV.12|) and (IIV.22|) . we obtain the magneti- 
zation ^A for the quasi-free electrons jTsj ] 



M = 2n 



a4 



knT 



m 



cos(27r£F /'huJc) 
^ smh{27r^ MkBT/heB) 



(V.40) 



where /iB is the Bohr magneton. The magnetic suscepti- 
bility X is defined by the ratio 



X 



M 
B ' 



(V.41) 



VI. 



TRANSPORT DIAMOND AND ZERO 
CURRENT ANOMALY 



We are now ready to discuss the TD and ZCA observed 
by Studenikin et al. [l|. Ref. Q, Fig. 1 is reproduced in 
Fig. [3l The outstanding features are: 

(A) The differential resistance r = dV/dl exhibit the 
SdH oscillations for higher DC, /dc — 50 /lA. The 
envelope of the SdH oscillations become smaller for 
weaker magnetic fields. 

(B) The background differential resistance for the SdH 



(C) The flat minima present at /dc = indicate a 
QHE. The flat minimum means a zero resistance 



R = V/I = 0. 



(VI.l) 
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FIG. 3: The differential resistance r of an InGaAs/InP Hall 
bar (width = 100 ^im) at different DC values, T = 270 mK. 
All curves except at Jdc ~ are shifted vertically by 0.25 kil 
for clarity. 



(D) The SdH maxima and the QHE minima both 
have the right-left symmetry with varying magnetic 
fields. 

(E) As the DC increases, the SdH maxima progressively 
become the QHE minima. 

Our interpretation is as follows. 

(A) The SdH oscillations are described by formula 
(|V.27|) . The oscillations are sinusoidal: 



cos(27r£F/?i'^c) ~ cos{2TTm* e-p / fieB) , 

and the envelope is represented by 

Trk^T 1 
£f smh{2TT^ M*kBT/neB)' 



(VI.2) 



(VI.3) 



The cyclotron mass m* appears in Eq. (jVI.2[) 
and the magnetotransport mass M* enters in 
Eq. (jVI.3|) . The two masses (m*,M*) correspond 
to the cyclotron motion and the guiding center mo- 
tion, respectively. We avoid the use of a Dingle 
temperature [Bj. 

(B) The background resistance (i?) averaged over the 
field B is zero: 



{R) 



0. 



(VI.4) 



This behaviour is in agreement with formula 
(jV.26p . It arises from the fact that there is no 
Landau-like term proportional to the squared mag- 
netic field B^ in the statistical weight W in 2D, see 
Eq. (|IV.24p . (There is no Landau diamagnetism in 
2D in contrast to the 3D case.) 

(C) The flat minimum meaning zero resistance R — 0, 
indicates the existence of a superconducting state. 
The superconducting state is stable with an energy 
gap. The supercurrents run with no scatterings 



by impurities and phonons. As is well known, the 
magnetic field is detrimental to the superconduct- 
ing state. If the excess magnetic field B* relative 
to the center field of the horizontal stretch exceeds 
a critical field, then the superconductivity is de- 
stroyed. The microscopic origin of this effect was 
explained in section 2. Briefly, the supercurrent is 
composed of the positively and negatively charged 
pairon-currents. The excess magnetic field B* gen- 
erates oppositely directed forces and breaks up ± 
pairons (Cooper pairs). 

(D) The magnetic field energy is quadratic in the ex- 
cess field B* , see Eq. (jll.lQp . which explains the 
right-left symmetry of the dcstroyment of the su- 
perconducting state. 

(E) The integer QHE occurring at the LL occupation 
numbers v — P — 1,2, . . . , have the quantized mag- 
netic fiuxes: 



BA 



1 
P 



(VI.5) 



where N^-^ is the fiuxon number at = P = 1. 
Hence the QHE has maxima at 



^ (VI.6) 



Equations (|V.27|) and (|V.29|) indicate that the re- 
sistivity r has minima when 



cos(27r£F/?ia;c) — 1, 

whose solutions are 
27reF 



huJr 



27rQ, g-1,2, 



We use ujc = eB/m*, ep = p|/2m*, 2ttpI 
and solve Eq. (jVI.8|) for B and obtain 



h 1 
B = -— ne 
e Q 



(VI.7) 



(VI.8) 



(VI.9) 



where n^^ is the electron density a,i v = 1. From 
Eq. pi. 171) . we obtain 



(VI.IO) 



Both P and Q are positive integers. Hence we find 
from Eqs. (lyTej) and (jVLOl) that the integer QHE 
maxima and the SdH minima occur precisely at 
the same magnetic fields B. Thus, the QHE min- 
ima progressively turn into the SdH maxima with 
increasing DC. 

Fig. 2, Ref. [l[ is reproduced in Fig.|4l The differential 
resistance r = dV/dl is plotted versus magnetic field (T) 
and direct current (/i A) . Diamond-shaped regions near 
SdH minima are called transport diamonds (TD). 
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FIG. 4: The differential resistance r of an InGaAs/InP Hall 
bar {w — 100 /^m) is plotted versus magnetic field and DC. 
The ZCA position is indicated by arrow. 
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FIG. 5: (a) The differential resistance of an InGaAs/InP Hall 
bar {w = 100 pim), T — 300 mK, is plotted versus magnetic 
field and DC; (b) the ZCA at two magnetic fields (1.58 T, 
1.70 T) indicated by vertical dashed lines in (a) for different 
temperatures. 



Our interpretation of the TD is a break-down of the 
superconducting QH state due to the excess magnetic 
field and the direct current. The direct current by itself 
generates a magnetic field, which is detrimental to the 
superconducting state. 

Fig. [5] is reproduced after Ref. [ij, Fig. 3. In (a) trans- 
port diamonds are shown, which are similar to those in 
Fig. [21 The temperature-dependence of the differential 
resistance is shown in (b) in the range (0. 253-1. 2K). 



The sharp dip in r vs. DC near DC = observed in 
Fig. [3] and Fig. [5ja) is called the ZCA, the narrow hori- 
zontal line indicated by an arrow in Fig. [31 Its tempera- 
ture behavior at i? = 1.58 T and B = 1.70 T is shown in 
Fig. [lib). 

The original authors [l| suspect that the origin of the 
ZCA arises from the Coulomb gap in the one particle 
density of states of interacting electrons. We propose a 
differing interpretation: 

Let us consider the case of the SdH at B — 1.58 T, 
the top figure in Fig. [S{b). At the low temperatures 
T = (0.253-1.2) K, the optical phonon population given 
by the Plank distribution function can be approximated 
by the Boltzmann distribution function: 



1 



(VI.ll) 



where £o is the longitudinal optical phonon energy (as- 
sumed constant) and fee the Boltzmann constant. The 
temperature dependence is exponential. The phonon 
population riph is rapidly changing with temperature 
and dominates. The resistance R is proportional to the 
electron-phonon scattering rate 7ph: 



R = a-^ 



7ph 



A, 



(VI.12) 



where Vph is the phonon speed, A the electron-phonon 
scattering cross section, and the phonon population 
given in Eq. (|VLlip . 

Studcnikin et al. 1] observed that the temperature de- 
pendence of the ZCA follows the Arrhenius law: 



7 oc e 

with the activation energy 



EA/kn = 1.3 K. 



(VI.13) 



(VI. 14) 



This value may correspond to the optical phonon energy 



(VI. 15) 



This finding supports our view that the temperature de- 
pendence of the ZCA arises from the electron-phonon 
scattering. 

We next consider the ZCA for the QHE at B = 1.70 T. 
This ZCA is also temperature-dependent. As the tem- 
perature decreases from 1.2 K to 0.253 K, the negative 
peak decreases in magnitude and its width becomes nar- 
rower. In the QHE under radiation, a supercurrent due 
to moving pairons condensed run in the upper (excited) 
channel and a normal current due to electrons run in the 
base channel. The resistance of the normal current is 
proportional to the electron-phonon scattering rate 7ph, 
as shown in Eq. (jVI.6p . Then, the phonon population 
approximately decreases exponentially at low tempera- 
tures (below 1.2 K). Thus the resistance decreases as the 
temperature T is lowered, which explains the observed 
temperature dependence. 
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The QHE at zero DC is destroyed either by increas- 
ing excess magnetic fields or by increasing DC-induced 
magnetic fields. But the ZCA indicates the destroyment 
is sharper for the case of increasing DC. This difference 
should arise from the direction of the magnetic field. The 
DC running along the sample length is likely to be inho- 
mogeneous, stronger at the outer edge. Then the super- 
conductivity is destroyed at the edges first according to 
Silsbeeb rule. On the other hand, the applied magnetic 
field alone should keep the current homogeneous. 

Studenikin et al. [l| observed essentially same TD and 
ZCA in hetero junction GaAs/AlGaAs. In particular, the 
QHE minima progressively turn to the SdH maxima as 
DC increases, and the ZCA is sharp near DC ~ 0. The 
same theory applies here. The authors thank Dr. S. Stu- 
denikin for enlightening discussions. 



Appendix A: DERIVA TION OF EQS. (flLS]) AND 

Dropping the "holes" from the Hamiltonian H in 
Eq. (|II.2|) . we obtain 



q 



k q 

- Vo 



k' 



where we suppressed the "electron" and spin indices. Us- 
ing the anticommutation rules pi.6p . we obtain 



■^C' ^kq - (^lk-hq/21 + '^|!.k+q/2|) ^kq 

- ^0 E' ^I'q (l " "k+q/2 - ^-k+q/s) ' 

(A.2) 



k' 



The Hamiltonian %c is bilinear in (_B, B"^), and can there- 
fore be diagonalized exactly: 



(A.3) 



where is the energy and cf)^ the annihilation operator. 
We multiply Eq. (jA.2p by (j)^ from the right, take a grand 
canonical ensemble average, denoted by angular brackets, 
and get 

w;^*^(k,q) = (e|k-Hq/2| +£|!.k+q/2|) *M(k,q) 



1 - /f (e|k'-|-q/2|) - /f (^l^-k' 



-Hq/2| 



(A.4) 



where (np) = fpi^p) is the Fermi distribution function. 
The reduced wavefunction 



(A.5) 



can be regarded as the mixed representation of 
the reduced density operator n defined through 

(k',q'|n|k,q) = ^^k^q^k'^q'^ The fc-boson energy 

can be specified by {Nj^, q), and it will be denoted by Wq 
since it is A^L-independent. As T — ^ 0, fpi^p) 0. Drop- 
ping the fluxon energy and replacing q/2 by q, we obtain 
Eq. (jlLSp . We solve this equation, assuming ep ^ Tic^D- 
Using a Taylor series expansion, we obtain Eq. (jll.6p to 
the linear in q. 



Appendix B: DERIVATION OF EQ. (flLT)) 

The BEC occurs when the chemical potential /i van- 
ishes at a finite T. The critical temperature Tc can be 
determined from 

n = {27:hr^ J d^p [e^=^ - l] (3, = {k^T,)-\ 

(B.l) 

After expanding the integrand in powers of e and 
using e = cp, we obtain 



n = 1.654(27r)-l(fcB^c/?^c)^ 
yielding Eq. (III.7I) . 



(B.2) 



Appendix C: STATISTICAL WEIGHT FOR THE 
LANDAU STATES 

The statistical weight W for the Landau states in 2D 
will be calculated in this appendix. We write the sum in 
Eq. pLiel) as 

OO 

2 ^ e(e - (2n + l)7r) = e(e - tt) + V(e; 0), (C.l) 

n=0 



'>p{e]x)= ^ e(e-7r-27r|7i + x|). (C.2) 

n— — OO 

Note that is periodic in x and can therefore be 

expanded in a Fourier series. After the Fourier expansion, 
we set X = and obtain Eq. (IC.ip . By taking the real 
part (Re) of Eq. (|C.1|) and using Eq. (jIV.20|) . we obtain 



Re{Equation ((Cl|) } = - / dre(e - r) 

Jo 

+ -y"(-lY dTQ(e-T) cos lyr, (C.3) 

where we assumed e = 2Tr E/TiUc 3> 1 and neglected tt 

against e. The integral in the first term in Eq. (jC.3P 
yields e. The integral in the second term is 

I dTQ{e — t) cosvT = — svuvt. (C.4) 
Jo ^ 
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We then obtain 



Re{Equation (ICAI) } -e + - V '— smi^e. (C.5) 



Using Eqs. pV.20j) and (|C.5j) . we obtain 



2 ^ (-1)'' . (2ttvE\ 



Chujc {-\ + Chujc- — sin 



(C.6) 



Appendix D: DERIVATION OF EQ. (|V.21|I 

Let us consider an integral on the real axis 

dx , z — X + iy and a, i? > 0. 

-R + 1 

(D.l) 

We add an integral over a semicircle of radius R in the 
upper z-plane to form an integral over a closed contour. 



We then take the limit: i? — >■ oo. The integral over the 
semicircle vanishes in this limit if a > 0. The integral on 
the real axis, I{y , a, oo), becomes the desired integral in 
Eq. (jV.21l) . The integral over the closed contour can be 
evaluated by using the residue theorem. Note that (e^ + 
1)^^ has simple poles at z = 7ri, 37ri, . . . , (2n — l)7ri, .... 
We may use the following formula valid for a simple pole 



at z — Zj '. 



Res{p(z)/(j(z), Zj} = p{zj)/q'{zj), 



(D.2) 



where p{z) is analytic at z — zj, and the symbol Res 
means a residue. We then obtain 



I {a, oo) = 27ri • ^ Res < z ' ~ ~ 



n=l 



°° gia[(2n~l)7ri] 

27ri • > — - — ... . — -27ri— 

/ ^ g(2n — l)in ^ 
n—1 

TT 1 



1 sinhaTT 



(D.3) 
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